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Abstract 



We report high-resolution Compton profiles (CP's) of Al along the three 
principal symmetry directions at a photon energy of 59.38 keV, together with 
corresponding highly accurate theoretical profiles obtained within the local- 
density approximation (LDA) based band-theory framework. A good accord 
between theory and experiment is found with respect to the overall shapes of 
the CP's, their first and second derivatives, as well as the anisotropies in the 
CP's defined as differences between pairs of various CP's. There are however 
discrepancies in that, in comparison to the LDA predictions, the measured 
profiles are lower at low momenta, show a Fermi cutoff which is broader, and 
display a tail which is higher at momenta above the Fermi momentum. A 
number of simple model calculations are carried out in order to gain insight 
into the nature of the underlying 3D momentum density in Al, and the role of 
the Fermi surface in inducing fine structure in the CP's. The present results 
when compared with those on Li show clearly that the size of discrepancies 
between theoretical and experimental CP's is markedly smaller in Al than 
in Li. This indicates that, with increasing electron density, the conventional 
picture of the electron gas becomes more representative of the momentum 
density and that shortcomings of the LDA framework in describing the elec- 
tron correlation effects become less important. 
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I. INTRODUCTION 



In a Compton scattering experiment one measures the so-called Compton profile (CP), 

J{Pz) = J J p{v>)dpxdpy, (1) 

where p(p) is the ground state electron momentum density. In an independent particle 
model the momentum density is given by 

p(p) = (27r)-3 Y, \J^{r) exp(«p • r)dr\^, (2) 



where ip{r) denotes the electron wave function.ETQ The summation in Eq. (|^) extends over all 
occupied states. The Compton profile, J{j>z), thus contains signatures of the Fermi surface 
breaks and correlation effects in the underlying three dimensional momentum distribution 
p(p). Since Fermi momenta pf are typically ~ 1 a.u., a high momentum resolution of ~0.1 
a.u. is essential in the experiment for delineating Fermi surface related fine structure in the 
CP. 

High-resolution Compton studies have recently been reported on Li,BQ Be,Era V,li^ and 
Cu.EJ In all these cases, careful comparisons of the shapes of the absolute valence elec- 
tron CP's, and the structure in the first and second derivatives as well as the directional 
anisotropics of the CP's, have been made with the corresponding parameter-free theoretical 
predictions based on the use of the LDA. A similar investigation of Li-rich LiMg disordered 
alloys where disorder effects were treated using the mean field KKR-CPA approach has also 
been carried out.li3 In this way, the band-theory based LDA approach has been shown to 
provide a remarkably accurate description of many aspects of the momentum density associ- 
ated with the quantum mechanical electronic ground state, including the characteristic fine 
structure induced by the Fermi surface. More exciting however is the fact that the afore- 
mentioned comparisons have for the first time clearly established the presence of systematic 
deviations between theoretical and experimental momentum densities. In Li, the experi- 
mental break Z^ in the momentum density at pf appears to be very small, nearly zeroQ; 



if so, this is very far from the results of electron gas calculations stretching over the last 
several decades.E3 In Be, the latest Compton dataH indicates anisotropic electron correlation 
effects outside the scope of much of the existing theoretical work which is based on treating 
properties of the homogeneous electron gas-li^rOJ For these reasons, a renewed interest in the 
problem of correlation effects on the momentum density beyond the LDA is natural,cHnH3 
although much further work is necessary for developing an approach of wide applicability in 
metals and alloys. 

Bearing these considerations in mind, there is strong motivation for undertaking high- 
resolution Compton studies of other systems. Our choice of Al in this connection is an 
especially good one because Al is trivalent and, therefore, it extends the range of electron 
densities investigated so far via high-resolution Compton. Correlation effects are of course 
expected to become less important with increasing electron density as the kinetic energy 
dominates. Also, Al has been the traditional touchstone of a free-electron-like metal with 
a nearly spherical Fermi surface (viewed in the extended zone). Neither a high-resolution, 
high-statistics Compton measurement, nor a band theory computation of high accuracy in 
order to identify Fermi surface related fine structure in the CP's of Al is currently available 
in the literature. E3'E3 The goal of the present work is to fill this gap and determine the 
extent to which the LDA describes the momentum density in Al. The existing Compton 
data on Al consists essentially of a number of measured CP's using 7-ray sources and solid 
state detectors at low momentum resolution.BZrEl Quite some time ago, Shiotani et a/.EiJ 
obtained the [111] CP of Al at a momentum resolution of 0.08 au, but did not investigate 
the anisotropy of the CP or the Fermi surface signatures therein. 

An outline of this article is as follows. In the next section we describe the experimental 
procedures. Section III gives pertinent details of computations. In section IV the experi- 
mental CP's are analyzed in the light of the band theory predictions as well as a number 
of other model computations. The Compton results are also compared briefly with closely 
related positron annihilation spectra. Section V summarizes our main conclusions. 



II. EXPERIMENT 

Single crystals of Al with surface normals oriented along the [100], [110] and [HI] direc- 
tions were used. The thickness of the crystals was about 2 mm. The reader is referred to 
Sakurai et alrn for details of our Compton spectrometer, and to Tanaka et alx^ and Itou et 
a/.B for our data processing procedures. Briefly, the spectrometer consists of a Cauchois-type 
bent-crystal analyzer of Si(422) with an image plate serving as a position sensitive detector. 
The scattering angle is 160°. The synchrotron radiation X-rays from a multipole wiggler 
installed in the 6.5 GeV Accumulation Ring at the National Laboratory for High Energy 
Physics are monochromatized by a quasi-doubly bent monochromator to 59.38 keV with an 
energy resolution of about 80 eV. The overall momentum resolution is estimated to be 0.12 
a.u. The double Compton scattering events were simulated via the Monte Carlo program of 
Sakaia; the integrated intensity of the double scattering events was found to be 10% of the 
single scattering events. 

The statistical error of each datum point, given by cr = V^+ 0.003A^, is estimated to be 
less than 0.3 %. Since the data points are not measured equidistantly they are interpolated 
onto an equidistant mesh of 0.02 a.u. using simple linear interpolation. The data were 
numerically differentiated according to y'{i) = \{[y{i + 1) — y{'i)]/[x{'i + 1) ~ ^^(i)] + [y{i) — 
y{i — l)]/[x{i) —x{i — 1)]}. No further smoothing or filtering was applied. The interpolation 
and differentiation cause some statistical correlation between the data points. 



III. COMPUTATIONS 

The band structure problem was solved within the all-electron charge self-consistent KKR 
framework without any free parameter. Exchange-correlation effects were incorporated using 
the von Barth-Hedin local spin density (LSD) approximation.l2a The lattice constant was 
computed to be 7.6534 a.u. by minimizing the total energy; for comparison, the experimental 



lattice constant at room temperature is 7.6559 a.u. The self-consistent crystal potential 
was obtained by iterating the KKR cycles using an elliptic contour with 48 points in the 
complex energy plane. The final charge density is self-consistent to an accuracy of about 
10~^ electrons and the Fermi energy to 10^^ Ry. An angular momentum cutoff Imax = 2 was 
employed. A free-electron-like Fermi surface was found with Fermi radii fcioo = 0.9246, knQ = 
0.9255, and km = 0.9292 a.u.; the free-electron value would be 0.9254 a.u. 

The CP's were obtained by first evaluating the three-dimensional momentum density 
p(p) in terms of the momentum matrix element of the KKR Green's functionEjE^ over a 
fine mesh of 48 x 4851 x 177 p-points, covering momenta up to Pmax ~ 5 a.u. This mesh 
involves 4851 k-points in the 1/48-th irreducible part of the Brillouin zone with each k-point 
translated into 177 p-points by adding reciprocal lattice vectors; the factor of 48 takes into 
account the symmetry operations of the cubic point group. The CP's can then be com- 
puted accurately by evaluating the two-dimensional integral of Eq. (|I]) using a generalized 
linear tetrahedron method.Ey The final CP's have been calculated over a momentum mesh 
containing 151 Pz points in the range 0-3 a.u. along each of the three measured directions. 
The accuracy of the computed profiles is about 1 part in 10^. A similar integration tech- 
nique has been used in our earlier studies of high- resolution CP's of various metals and 
alloys.3ErBllllE3El'E3. The Lam-Platzman correction^ to the CP's was computed using the 
occupation number density of the uniform electron gas. 



IV. RESULTS AND DISCUSSION 

Figure |l] shows the measured and computed CP's of the valence electrons along the 
[100], [110] and [111] directions; the theoretical CP's are convoluted with a gaussian which 
represents the experimental resolution of 0.12 a.u. FWHM. The experimental valence CP's 
have been obtained by subtracting the theoretical core CP's from the measured profiles. In 
this connection, we used the solid-state core wavefunctions which refiect the shght overlap 



of the 2p core states in Al. The first and second derivatives of the valence profiles have been 
obtained by numerical differentiation. 

In examining the overall shape of the CP's in Fig. ^ one notes that the experimental 
points are lower at low momenta compared to the calculated values. We emphasize that this 
does not imply that the measured 3D momentum density is lower than the theoretical one 
at all momenta. To see this, recall thatch 

1 (PJaviv) 



P(0) 



(3) 

p=0 



27r rfp2 

where Jav. (p) denotes the directionally averaged CP which in a cubic crystal may be reason- 
ably approximated byn3 

JavXp) = (l/35)[10Jioo(p) + 16Jiio(p) + 9Jiii(p)]. (4) 

The bottom row in Fig. |l| shows that the differences between the experimental and theoretical 
second derivatives at p^ = are well within the error bars. In view of Eq. (^, this indicates 
that the underlying 3D distributions are not significantly different at Pz = 0. In fact, this 
result implies that the measured momentum density must be smaller than the theoretical 
one at momenta approaching the Fermi momentum pf. This is also borne out by the first 
derivatives shown in the central row of Fig. |l] which begin to show differences between 
experiment and theory only above Pz ~ 0.3 a.u. 

Further insight is provided by Fig. ^ which shows the spherically averaged 3D momentum 
density defined by 

P^^^Xp) = -il/2np){dJav./dp), (5) 

where Jav.(p) is obtained via Eq. (^. The oscillations in Pav.(p) at small momenta reflect 
partly the large (correlated) error bars due to the division of the small derivative by small 
values of p, and partly the (spherically averaged) effect of Brillouin zone-face interactions to 
be discussed below. In any event. Fig. |^ makes it clear that the experimental 3D momentum 
density lies below the theoretical predictions as one approaches pj, and that the situation 
reverses itself above pf. 
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The 2nd derivatives J" in Fig. [| all show a peak at pf. There is good agreement between 
theory and experiment as to the position of the peaks, i.e., the value of pj, but the measured 
peaks are all lower and broader than the theoretical predictions. Although the shapes of 
the J" peaks reflect the complex interplay between the effects of experimental resolution, 
electron correlations and lattice potential on the Fermi cut-off in the momentum density, it 
is evident from Figs. || and § that the measured distribution possesses a tail higher than 
the theory beyond pf. 

These discrepancies between theory and experiment are similar to those reported earlier 
in LiBQ and other metalglrEHllI and have their origin in the electron correlation effects beyond 
the LDA which are not treated properly within our theoretical framework. Such correlations 
are expected to cause (relative to the independent particle model) a decrease of the momen- 
tum density as one approaches p/, and a tail at momenta greater than pj; as indicated above, 
both features are qualitatively visible in our comparison between theory and experiment. 
Notably, the deviations from LDA theory are smaller in Al than in Li. For example, the dif- 
ference between the theoretical and experimental valence profiles at p^ = is approximately 
16 % for Li and 4.5 % for A1,e3 and the width of the peak at pf in the second derivatives is 
0.23 a.u. in Li and 0.15 a.u. in Al; thus, the "blurring" of the Fermi cutoff is more severe 
in Li than in Al. These characteristic differences between Li and Al are partly related to 
the difference in the electron density of the two metals. The electron density in terms of r^ 
(the standard parameter for the volume per electron of valence electrons), is 3.21 for Li and 
2.12 for Al. Therefore, the bare Coulomb interaction is more effectively screened in Al than 
in Li. As shown by a variety of treatments of the homogenous interacting electron gas, as 
the electron density increases, the kinetic energy dominates, and the momentum density is 
described more closely by the free-electron rectangular distribution with a step-wise cutoff 



at p^.oilHiZhlHl 

Figure 3 considers the effect of the isotropic Lam-Platzman (LP) correction on the [111] 
CP; results along other directions are similar and are not shown in the interest of brevity. 
The theoretical curves in Fig. |^ include the LP-correction, while those in Fig. do not. A 
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comparison of Fig. |^ with the last column of Fig. [l| shows that, although the inclusion of 
the LP-correction improves things, much of the discrepancy between theory and experiment 



still remains. Interestingly, Ref. ^ has recently analyzed the correlation correction to the 
CP's of Al in terms of a model which involves the break Z^ in the momentum density at pf 
as the only free parameter. By adjusting Z^, Ref. |2^ finds that the discrepancy between the 
LDA predictions and the measurements can be essentially removed for a Zk value between 
0.7 and 0.8, in reasonable accord with the corresponding theoretical values from various 



authors which are scattered between 0.76 and O.SS.EjIiil There is no inconsistency between 



the present results and those of Ref. ^. To see this relationship, recall that the standard 
LP-correction is defined viaH 

Apip) = J dhp{v) [p^^^b, r,(r)) - p^^(p, r,(r))] , (6) 

where the integral extends over the Wigner-Seitz cell. The expression within the square 
brackets gives the difference between the momentum densities of the interacting and non- 
interacting homogeneous election gas (denoted by the superscripts 'INT' and 'NI') evaluated 
at the local density p(r) of the physical system and rs(r) is the corresponding electron density 
parameter. Equation (P) thus attempts to take into account inhomogeneities in the electron 
gas, whereas the semi-empirical model of Ref. ^ replaces the integrand by its value at the 
average electron density in Al. The matter is quite subtle, and further work is necessary in 
order to develop a satisfactory treatment of correlation effects on the momentum density in 
solids. 

If the momentum density within the Fermi sphere were flat and smooth, the first deriva- 
tive dJ{pz)/dpz shown in the middle row of Fig. |1| would be a straight line up to the cutoff 
at the Fermi momentum. However, at momenta less than the Fermi radius some structure 
is visible notably in the [111] and the [100] derivatives. In this connection, we note that the 
Fermi sphere overlaps with Umklapp Fermi spheres centered on the (111) reciprocal lattice 
points around the W-points in the Brillouin zone. For example, the hexagonal zone face 
contains six W-points which all project at (|, |, 1), i.e., the point pz = 0.71 a.u. on the [111] 
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axis. Similarly, four W points in the first Brillouin zone project at 0.41 a.u. and another four 
at 0.82 a.u. on the [100] direction. Interestingly, the experimental as well as the theoretical 
derivatives contain structure around 0.7 a.u. in the [111] and 0.4 a.u. in the [100] CP. This 
indicates the importance of the k states near the W points with respect to the fine structure 
in the Al CP's. Incidentally, a structure similar to the wiggle around 0.4 a.u. in the [100] 
derivative has been observed in positron annihilation ID-ACAR measurements by Okada 
et a/.E3 and 2D-ACAR measurements by Mader et alM who also ascribed it to zone-face 
interactions around the W-points. 

The directional differences, shown in Fig. ^, are a measure of the anisotropy. Although 
the maximum difference is about 1 % of the peak value of the profile itself, they show 
definite structures which can have several origins. Firstly, for different crystal orientations 
the plane of integration in Eq. (|l]) sweeps differently through the Umklapp Fermi spheres 
centered at the reciprocal lattice points in the higher Brillouin zones. Secondly, the Fermi 
surface is slightly distorted from a sphere, as witnessed by the different Fermi radii given 
above, while, thirdly, band-structure effects such as a p dependence of the momentum density 
within the Fermi spheres and interactions of the electron bands with the Brillouin zone faces 
with consequent distortion of the wavefunctions will also contribute to the anisotropy. The 
importance of the first point can be readily studied using a simple model of a spherical free- 
electron-like Fermi surface surrounded by seven shells of similar Umklapp Fermi surfaces. 
The momentum density within each Fermi surface is assumed to be constant and given by 
the square of the corresponding Fourier component of the electron wave function at TiB The 
CP for a given direction then consists of a superposition of inverted parabolas, centered at 
the projections on that direction of the reciprocal lattice points. The height of each parabola 
is proportional to the momentum density within the corresponding Fermi sphere, while its 
cut-off points are found by adding or subtracting pf from the projected center. Figure |^ 
shows the directional differences thus obtained. The positions of the cut-off points have 
been indicated by the arrows at the bottom of the graph, together with a symbol which 
denotes the direction of projection (a = [100], b = [110], c = [111]) and the coordinates 
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of the center of the Fermi sphere. It should be noted that many Umklapp Fermi spheres 
coincide in projection and therefore these coordinates are not unique; the simplest set has 
been noted. The analysis in Fig. ^ shows that much of the important structure in the 
directional differences stems from the < 111 > Umklapp contributions; the other Umklapps 
play a less important role. 

A comparison of Fig. ^ with the calculated differences in Fig. ^ shows an overall quali- 
tative correspondence in the succession of positive and negative peaks. On a more detailed 
scale, however, there are significant differences which have their origin in the other factors 
mentioned above. Notable examples are the peaks around 0.85 a.u. in the J[iii] — J[ioo] 
and J[iii] — J[iio] directional differences in Fig. ^ which have no clear counterpart in Fig. |^. 
Kubo et alx3 have ascribed these features to the fact that in the [111] direction the actual 
Fermi surface bulges out beyond the free-electron Fermi sphere in the second Brillouin zone 
while there is a contraction in the third zone. This will strongly affect the [111] profile but 
not so much the other two. Our simple free-electron model of course does not contain this 
Fermi surface distortion effect. The calculated curves in Fig. ^, on the other hand, include 
all of these factors and reproduce the essential characteristics of the measured differences, 
although some discrepancies remain. It may be noted that non-locality of the exchange and 
correlation potential in Li reduces the Fermi surface anisotropyHi^ and thus would affect 
the anisotropy of the CP's. In this vein, lattice vibrations would reduce the Umklapp con- 
tributions and hence the anisotropy of the momentum density. In how far such effects can 
explain the residual discrepancies in Fig. |^ remains unclear. 

In principle, the anisotropy in the momentum density may be obtained approximately by 
expanding both the momentum density and the CP's into lattice harmonics, and establishing 
the relation between the expansion coefficients for the momentum density and those for the 
CP's.E3 Actually, Eq. (^) represents the / = term in such a scheme. However, we have not 
attempted to analyze our data along these lines since the number of measured profiles is not 
large enough. 

Additional information may be gained from a comparison of CP's with the correspond- 
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ing results of positron annihilation measurements. Both experiments probe the momentum 
density - in positron annihilation one measures the momentum density of the annihilating 
electron-positron pair whereas in Compton only the electron momentum density is involved. 
In Fig. ^ the first derivative of the one-dimensional angular correlation of positron annihi- 
lation radiation (ID- AGAR) profile measured by Okada et a/.E3 for the [111] orientation is 
compared with the corresponding CP. The momentum resolution of the ID-ACAR is 0.11 
a.u. which is almost the same as that of the present CP's. Since the ID-ACAR was area 
normalized to the theoretical ID-ACAR calculated by Kubo et al.p^ the peak height at 
P2 = is almost the same as that of the present CP. The slope at the Fermi momentum is 
steeper in the ID-ACAR than in the CP, which is direct evidence for enhancement of the 
annihilation of positrons with the s — p electrons near the Fermi energy predicted first by 
KahanaE3 on the basis of an interacting electron gas model. Also, the correlation tail for 
p > Pf in the ID-ACAR is weaker than its counterpart in the CP as a result of the partial 
cancellation of electron-electron and positron-electron correlation effects.^ Finally, the fine 
structure at 0.2 a.u. and 0.7 a.u. is more pronounced in the ID-ACAR than in the CP. This 
points to less correlation- induced smearing in positron annihilation compared to Compton 
scattering. 



V. SUMMARY AND CONCLUSIONS 

We have measured the Compton profiles (CP's) of Al along [100], [110] and [111] direc- 
tions at a photon energy of 59.38 keV and a momentum resolution of 0.12 au. Parallel, highly 
accurate all-electron computations have been carried out within the LDA-based band-theory 
framework. Comparisons between theory and experiment at the level of the shapes of the 
CP's, structure in the 1st and 2nd derivatives of the CP's, and anisotropics obtained by 
taking differences between three pairs of CP's, all show a good level of accord. However, 
there are discrepancies as well. In comparison to the LDA predictions, the measured profiles 



13 



are lower at low momenta, show a Fermi cutoff which is broader, and display a tail which 
is higher at momenta above the Fermi momentum. The inclusion of correlation effects in 
the LDA via the standard isotropic Lam-Platzman correction improves things slightly, but 
the essential discrepancies remain. A model analysis in terms of directionally averaged CP's 
allows us to get a handle on the 3D momentum density of Al; in this way, we adduce that 
the experimental 3D density near p=0 does not differ significantly from LDA predictions 
even though the CP's do. In this vein, CP's are computed using a model 3D distribution in 
which free electron spheres with appropriate weights are placed on reciprocal lattice points 
(extending to seven shells around a central sphere) to represent the higher momentum com- 
ponents in the electronic wave functions; the results show that a significant amount of fine 
structure in the CP's is induced by these higher momentum components and by k- states 
near the W-points in the Brillouin zone where the free electron spheres overlap. 

The present results when compared with those reported earlier on Li show clearly that 
the size of discrepancies between theoretical and experimental CP's is markedly smaller in 
Al than in Li; in particular, theoretical and experimental profiles at p^ = differ by 16% in 
Li but only by 4.5% in Al, and the peak width in the 2nd derivative at pf is 0.23 au in Li 
but 0.15 au in Al. It is thus clear that, with increasing electron density, the conventional 
picture of the electron gas becomes more representative of the momentum density and that 
shortcomings of the LDA framework in describing the electron correlation effects become 
less important. Finally, we compare briefly our [111] CP with the positron-annihilation (ID- 
ACAR) measurements of Okada et al., and show that in the case of positron- annihilation 
the Fermi cut-off is sharper and that there is less correlation induced smearing of structures 
in the ACAR spectrum. 
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FIGURES 
FIG. 1. Top: Measured and computed Compton profiles of Al along the [100], [110] and [111] 

directions. Theoretical profiles (solid lines) have been broadened to reflect experimental resolution. 

Middle: First derivatives of the measured and computed profiles. Bottom: Second derivatives of 

the measured and the computed profiles. 

FIG. 2. Theoretical (solid curve) and experimental (dashed curve) directionally averaged 3D 
electron momentum density obtained via Eqs. ^ and (|5|). 

FIG. 3. Same as the last column of Fig. |l], except that here the theoretical curves in all cases 
include the Lam-Platzman correction. 

FIG. 4. Measured and computed directional difference profiles for three different pairs of 
directions. 

FIG. 5. Directional difference profiles calculated for a simple quasi- free-electron model of Al in 
which the CP is given by a superposition of parabolic contributions centered at various reciprocal 
lattice points (see text). The arrows at the bottom indicate the positions of the cut-off points (each 
parabola has two cut-off points; the other one lies outside the graph). The directions of projection 
are indicated by a (=[100]), b (=[110]), and c (=[111]), while the subscripts denote the coordinates 
of the reciprocal lattice points involved. 000 denotes the cut-off of the central Fermi surface, i.e., 
the Fermi radius pf. 

FIG. 6. First derivative of the ID- AGAR spectrum (open circles) along the [111] direction read 
off from Ref. ^ is compared with the derivative of the [111] Compton profile shown in Fig. |^. 
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